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In this work we have studied a new functional for the correlation energy obtained from the exact- 
exchange (EXX) approximation within time-dependent density functional theory (TDDFT). Cor- 
relation energies have been calculated for a number of different atoms showing excellent agreement 
with results from more sophisticated methods. These results loose little accuracy by approximat- 
ing the EXX kernel by its static value, a procedure which enormously simplifies the calculations. 
The correlation potential, obtained by taking the functional derivative with respect to the density, 
turns out to be remarkably accurate for all atoms studied. This potential has been used to calcu- 
late ionization potentials, static polarizabilities and van der Waals coefficients with results in close 
agreement with experiment. 
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I. INTRODUCTION 

Time-dependent density functional theory (TDDFT) 
provides a promising and rigorous framework for treat- 
ing interacting many-electron systems within a reason- 
able computational cost. Although the overall aim of 
TDDFT is to describe physical phenomena associated 
with excited states, strong connections between the latter 
and ground-state properties make it possible to obtain an 
improved description of the ground state through the use 
of TDDFT. 

In many interesting cases, it is sufficient to treat sys- 
tems under the influence of external perturbations weak 
enough to allow for a description in terms of linear 
response. Within this realm, the basic quantities of 
TDDFT are the ground-state exchange-correlation (XC) 
potential and the corresponding XC kernel (/xc)- 1 In pre- 
vious work we have studied these quantities at different 
levels of approximation. The XC potential has been cal- 
culated for atoms at the level of the GW approximation 2 
and the XC kernel has been presented within the exact- 
exchange (EXX) approximation. 3-5 In all these previous 
publications the necessary formulas have been derived 
from the variational formulation of many-body perturba- 
tion theory. 6 ' 7 The virtue of this approach is that the ob- 
tained results are guaranteed to obey many conservation 
laws and sum rules like, e.g., the /-sum rule, of impor- 
tance to the calculated optical spectra. In the present 
work we have deviated from this path and instead fol- 
lowed an approach originally suggested by Peuckert. 8 

Having some approximation for the density-density re- 
sponse function of a many-electron system one can obtain 
an expression for the total ground-state energy by using 
the Hellmann-Feyman theorem applied to the strength 
of the inter-particle Coulomb interaction. 9 ' 10 And, by 
means of TDDFT, from any approximation for the XC 
potential and the corresponding XC kernel one can obtain 
an approximation to the density response function of the 
system. The XC part of the resulting total energy can 



then be differentiated once with respect to the density 
to yield a new approximation for the XC potential and 
then twice to yield a new approximation for the XC ker- 
nel. From these results we can obtain a new total energy 
which again can be differentiated to obtain a new poten- 
tial and a new kernel, and so on. Of course, in practice 
the resulting expressions quickly become unmanageable 
and the proposed iterations have to be limited to one or 
two steps. 

In previous work we have seen that the total energy ob- 
tained via the Hellman-Feynman theorem applied to the 
EXX approximation is very accurate giving errors of the 
order of 5% in the resulting correlation energies. 3 These 
results inspired us to believe that a differentiation of the 
corresponding expression for the XC energy with respect 
to the density might give rise to a very accurate XC po- 
tential and XC kernel. And this is, indeed, what we have 
found in the present work, at least as far as the poten- 
tial is concerned and, then, within the approximations 
we have been forced to do in order to obtain tractable 
expressions. 

In the present context, the XC kernel of the EXX ap- 
proximation (EXXA) has the convenient property of be- 
ing linear in the Coulomb interaction allowing us to carry 
out the integration over the strength of the Coulomb in- 
teraction analytically. The result is a closed expression 
for the XC energy. This expression explicitly contains the 
EXX kernel / x , which give rise to numerical difficulties 
in the later process of differentiating the XC energy with 
respect to the density. Motivated by the fact that the 
EXX kernel for He is independent of the density we have 
neglected this density dependence for all atoms. The va- 
lidity of this approximation must, of course, be verified 
independently but, encouraged by our excellent results, 
we have decided to leave this test to a future publication. 

We have also employed an additional approximation 
which we believe to be of even less consequence. The 
fact that the EXX kernel is frequency dependent for all 
systems but for a two-electron one leads to much longer 
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computational times. We have, however, seen that total 
energies are very insensitive to this frequency dependence 
and we thus recommend to neglect it. The frequency 
dependence will, of course, also affect the calculation of 
the XC potential except in the case of He. We have here 
assumed that also the potentials of the heavier atoms are 
relatively insensitive to this frequency dependence but we 
leave also the verification of this assumption to a future 
publication. 

The XC potentials that we obtain from the first itera- 
tion of the Peuckert procedure starting from the EXXA 
turns out to be much better than any approximate XC 
potentials perviously obtained and they are actually very 
close to the exact ones 11 where these are known. For a 
number of spherical atoms, we have used these new po- 
tentials to calculate total energies, static polarizabilitics, 
and van der Waals coefficients using the expression for 
the XC kernel obtained within the EXXA. In all cases 
we find excellent agreement with experiment. We thus 
conclude that we now have an affordable and well de- 
fined way of obtaining accurate results for ground-state 
properties and low lying excitations of many-electron sys- 
tems. 



II. CORRELATION ENERGY FUNCTIONAL 

A standard expression for the correlation energy can 
be obtained by introducing a fictitious Hamiltonian H x 
with a scaled Coulomb interaction \v and a local multi- 
plicative potential which guarantees that the density is 
constant for every value of the scaling parameter A. At 
A = 1, H x coincides with the fully interacting Hamilto- 
nian and at A = with the one of the non-interacting 
Kohn-Sham (KS) system. Using the Hellman-Feynman 
theorem one can show 9 ' 10,12 
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where Xs is the non-interacting KS response func- 
tion and x A is the scaled density response function. 
We have also used the short hand notation Txfg = 
J drdr'f(r, r')g(r',r) for any two-point functions / and 
g. Within TDDFT the function x A reads 



X X =Xs+Xs [A« + / X A C ] X X - 



(2) 



The scaled XC kernel / xc is a functional of the ground- 
state density and is defined as the functional derivative 
of the scaled XC potential u xc with respect to the density 
n. 

The simplest approximation to Eq. (1) is the random 
phase approximation (RPA), obtained with / xc = 0. The 
RPA has the advantage of allowing for an analytical eval- 
uation of the A- integral in Eq. (1). The result is 
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In the language of Feynman diagrams Eq. (3) is equal to 
an infinite summation of ring-diagrams. 

The RPA correlation potential v c is obtained as the 
functional derivative of Eq. (3) with respect to the den- 
sity. If we let V signify the total KS potential, G s the 
non-interacting KS Green function, Xs = —iG s G s and 
the functional derivative is conveniently obtained via the 
chain rule 
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The result is the well-known linearized Sham-Schliiter 
(LSS) equation 6 - 13 

J 2)« c (l)d2 = J A(3, 2; 1)E C (2, 3)d2d3. (4) 

Here, we have used the notation (ri,ti) = 1 etc. and 
introduced A(3,2;l) = -iG g (3, 1)G 3 (1, 2). The correla- 
tion part of the self-energy S c in the RPA is given by 



^ = iv X ™ A vG. 



where 
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Thus, the S c in the RPA coincides with the GW self- 
energy but evaluated at KS Green functions. For the 
purpose of obtaining the RPA potential Eq. (4) and Eq. 
(6) have to be solved self-consistcntly and has only been 
done so far for atoms 2 and in bulk Si, LiF and solid 
Ar. 14 The RPA is also sometimes called the linearized 
time-dependent (TD) Hartree approximation since x RPA 
is obtained by allowing the electrons to respond only to 
the perturbing potential plus the induced Hartree poten- 
tial. The next level of approximation is obtained by also 
including exchange effects, which leads to the RPAE ap- 
proximation or the linearized TD Hartree Fock (TDHF) 
approximation. Within TDDFT the same level of ap- 
proximation corresponds to the TDEXX approximation, 
in which the HF potential is replaced by the local EXX 
potential v x . The latter potential is obtained from the 
TD exchange version of the LSS equation which amounts 
to replacing, in Eq. (4), v c by v x and S c by £ x = ivG a , 
i.e., the HF self-energy. A variation of that equation 
yields an equation for the EXX response kernel / x : 



J Xs (l,2)/ X (2,3)x s (3,4)rf2d3 
= j <5S X (2,3) 



-A(3,2;l)d2d3 



+ 



5V(4) 

J A(l,2;4)A(2,3)G s (3,l)d2d3 
+ J G s (l,2)A(2,3)A(3,l;4)d2d3, 



(7) 



where A(2, 3) = S x (2, 3) - v x (2)6(2, 3). A full analysis of 
this kernel has been performed recently. It is important 
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to observe that the kernel / x is an implicit functional of 
the density through the KS orbitals and their eigenvalues. 
From Eq. (7) we see that the potential v x is an ingredient 
in the kernel / x but, as discussed above, v x is also an 
implicit functional of the density. For future reference we 
here note that an EXX kernel obtained for any density 
and used in Eq. (2) will generate an interacting response 
function which obeys the important /-sum rule. This 
is due to the fact that the procedure above guarantees 
that the kernel does not blow up at large frequencies, 
something that we proved in a previous publication. 3 The 
fulfillment of the /-sum rule is crucial for the construction 
of new functionals for the correlation energy, one of which 
we will now derive. 

We start by making the observation that the EXX ker- 
nel is linear in the explicit dependence on the Coulomb 
potential. Therefore, the A-integration in Eq. (1) can, 
just as in the case of RPA, be carried out analytically 
yielding the result 

(8) 

We note that Tr^Xs contains a singularity proportional 
to the Coulomb potential at the origin times the number 
of particles. This singularity must be cancelled by the 
first term in Eq. (8) which will only occur if the kernel 
remains finite at large frequencies, i.e., obey the sum rule. 

Correlation energies obtained from this functional was 
recently presented 3 but in a non self-consistent fashion 
using the EXX density for the evaluation. The results 
were excellent for many closed-shell atoms. A diagram- 
matic description of Eq. (8), up to second order, is given 
in Fig. 1. Many higher order terms do not have a strict 
diagrammatic representation but can be considered to 
simulate the higher order exchange diagrams. 

In order to obtain the correlation potential from Eq. 
(8) we need to differentiate the correlation energy with 
respect to the density. We then encounter the problem of 
differentiating the EXX kernel with respect to the den- 
sity, which is a rather cumbersome task. Such deriva- 
tives formally amounts to three-point correlation func- 
tions 5/ x (l, 2)/<5n(3). For He, however, this derivative is 
zero since / x = — ^v. Therefore, we have here chosen to 
see how well we can do by neglecting the density varia- 
tion of the EXX kernel also for larger atoms. With this 
assumption the same procedure used for going from Eq. 

(3) to Eq. (4) now leads to an equation similar to Eq. 

(4) with E c replaced by 

X c = iv X [v + f x ]G s , (9) 

where \ is the interacting response function in the 
TDEXX approximation (Eq. 2 at A = 1 and / xc = / x ). 
The evaluation of the correlation energy form Eq. (8) is 
relatively time consuming due to the frequency depen- 
dence of the kernel / x . 3 The computational cost would 
be substantially reduced it the kernel could be kept at 
its static value (/ x (r, r',w) w / x (r, r',0)) without seri- 
ously affecting the correlation energies. In Table I the 




FIG. 1: With / xc = / x , the correlation energy functional has 
a diagrammatic expansion up to second order. Higher order 
terms are simulating higher order exchange diagrams. 



correlation energies evaluated using the fully frequency- 
dependent EXX kernel (TDEXX) are compared to those 
obtained by using the static (=adiabatic) kernel (AEXX) 
and the difference is seen to be relatively small. Notice 
that in the case of He the static kernel is the full kernel. 
These results have encouraged us to evaluate the corre- 
lation potentials also using the adiabatic approximation. 



III. CORRELATION ENERGIES AND 
POTENTIALS 

When evaluating the correlation energy for a system it 
is natural to use the self-consistent density for that par- 
ticular system. From now on, the self-consistent density 
corresponding to the new functional in Eq. (8) is referred 
to as the RPAX density. Due to the stationary property 
of the total energy we expect, however, that evaluating 
correlation energies at a slightly different density will give 
almost the same result. 15 This is, indeed, the case as 
can be seen in Table I, which also seems to demonstrate 
that the total energy is a minimum at the RPAX den- 
sity. Clearly, the EXX kernel gives a large improvement 
over the too large RPA values. It also improves the MP2 
results, which are here the self-consistent results given in 
Ref. 16. The latter approximation also follows from Eq. 
(1) if X A is replaced by x A ~ Xs + ^Xs [v + / x ] Xs, or if the 
logaritm in Eq. (8) is expanded to second order. The use 
of the adiabatic approximation is seen to be less severe 
yielding energies of the same quality as those obtained 
with the frequency-dependent kernel. As noticed before, 
the values in TDEXX, or in AEXX for that matter, are 
very accurate for these systems. 

The XC potential is interesting in its own right. The 
highest occupied eigenvalue of the KS system exactly cor- 
responds to the ionization potential 18 and the larger part 
of the particle-conserving excitation energies consists of 
KS eigenvalue differences. 19 In Tab. II we present ion- 
ization potentials produced by different KS potentials 
(RPAX, RPA, MP2 and EXX). As noticed previously, the 
RPA values improve over EXX and are also better than 
the MP2 values. Here, we see that the RPAX potential 
yields an even further improvement giving excellent ion- 
ization potentials for all atoms. In Fig. 2 we plot the same 
correlation potentials for He, Be, and Ne. The He RPAX 
potential almost coincide with the exact potential 11 and 
the RPAX potential for Be is much closer to the exact 
than any other approximation we have tried - especially 



4 



TABLE I: Correlation energies from a few different approx- 
imations. For a consistent comparison with CI results the 
correlation energy is here defined as the difference between 
the total energy and the Hartree-Fock energy. The fourth 
decimal is shown in parentheses in order to compare different 
approximations, (a.u.) 



He 
Be 
Ne 
Mg 
Ar 



TDEXX 
EXX 



AEXX 
EXX 



TDEXX 
RPAX 



AEXX 
RPAX 



0.044(5) 
0.102(0) 
0.388(9) 
0.445(2) 
0.727(8) 



0.044(5) 
0.101(7) 
0.377(1) 
0.437(2) 
0.710(6) 



0.044(6) 
0.103(3) 
0.390(3) 
0.446(6) 
0.728(7) 



0.044(6) 
0.102(8) 
0.377(8) 
0.438(4) 
0.711(2) 



RPA 
RPA 



0.083 
0.181 
0.596 
0.681 
1.091 



MP2 
MP 2 



0.047 
0.124 
0.480 
0.514 
0.844 



CP 



0.0420 
0.0943 
0.3905 
0.4383 
0.7225 



'From Ref. 16. 
c From Ref. 17. 



in the outer region. The Ne potential is also very accu- 
rate, yielding even a qualitative improvement by better 
describing the 2s-shell giving rise to an extra structure in 
the RPAX, and the exact, correlation potential. 



IV. STATIC POLARIZABILITIES AND VAN 
DER WAALS COEFFICIENTS 

The static polarizability is defined according to the 
formula 



o(0) = — I zx(t,t',lj — 0)z'drdr', 



(10) 



and the van der Waals coefficient, or Cg-coefficient, be- 
tween ion A and B is given by 



o.a(zo;)o b {ioS)doj, 



(11) 



where 04(10;) is the dynamic polarizability of ion A cal- 
culated at imaginary frequencies. In previous work these 
quantities were calculated using the response function of 
the TDEXX approximation. It was then a natural choice 
to use the corresponding self-consistent EXX density in 
the evaluation. We then found, not so surprisingly, that 
our results closely resembled those of the TDHF approx- 
imation and were, therefore, not overly impressive. We 



TABLE II: Ionization potentials obtained from the highest 
occupied KS eigenvalue of different KS potentials, (a.u.) 



Atom 


EXX 


MP2 


RPA 


RPAX 


Exp. 


He 


0.918 


0.893 


0.902 


0.904 


0.904 


Be 


0.309 


0.357 


0.354 


0.340 


0.343 


Ne 


0.851 


0.657 


0.796 


0.787 


0.792 


Mg 


0.253 


0.302 


0.297 


0.282 


0.281 


Ar 


0.591 


0.558 


0.590 


0.577 


0.579 



interpreted this partial failure to be a consequence of a 
rather poor description of the ground state. We found, 
however, also that the results are rather sensitive to the 
density used in the evaluation. In the present work we 
have instead evaluated the same EXX formula for the po- 
larizability using the correlated density produced by our 
new correlation functional. The van der Waals energy is 
a pure correlation effect and one can argue that it should 
be evaluated at a correlated density and not just the ex- 
change density like, e.g., the RPA or the RPAX density. 
This turns out to have rather drastic effect on the actual 
values moving them much closer to the more accurate 
values found in the literature as seen in Tab. Ill and IV. 
A clear improvement is found when correlated densities 
are used in the evaluation with a slight edge for our new 
RPAX density. It is also noticed that when calculating 
Cg-coefficients the adiabatic approximation is sufficient 
as also observed previously. 24 
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FIG. 2: Self-consistent correlation potentials for He, Be 
and Ne in different approximations. The MP2 potential for 
Be is evaluated at the EXX density due to a variational 
instability. 16 The exact potentials are those of Umrigar et 
al. 11 
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TABLE III: Static polarizabilities calculated from \ m the 
RPA and the TDEXX approximation. The latter has been 
evaluated using different potentials (EXX, RPA, RPAX, and 
Exact), (a.u.) 



JXC : 


TDEXX 








RPA 


Litt. 




EXX 


RPA 


RPAX 


Exact 


RPA 




He 


1.322 


1.351 


1.348 


1.349 


1.225 


1.38 a 


Ne 


2.372 


2.577 


2.613 


2.555 


2.424 


2.67 a 


Ar 


10.74 


10.69 


10.94 




9.839 


11.08 a 


Be 


45.64 


40.09 


41.04 


40.49 


28.99 


37.7 6 


Mg 


81.66 


70.37 


71.67 




51.56 


71.35 6 



"Prom Ref. 20 
''From Ref. 21 



TABLE IV: van der Waals or Cq coefficients calculated from 
X in the RPA, the AEXX and the TDEXX approximation. 
The latter has been evaluated using different potentials (EXX, 
RPA, RPAX, and Exact), (a.u.) 



/xc • 


TDEXX 






AEXX 


RPA 


Litt. 


«x C : 


EXX 


RPA 


RPAX 


Exact 


RPAX 


RPA 




He 


1.375 


1.414 


1.411 


1.411 


1.411 


1.206 


1.458° 


Ne 


5.506 


6.091 


6.191 


6.021 


6.161 


5.523 


6.383° 


Ar 


61.88 


61.27 


63.19 




63.11 


53.69 


64.3° 


Be 


282.8 


226.7 


235.3 


231.5 


236.5 


142.0 


214 6 


Mg 


767.5 


617.8 


634.3 




632.3 


385.6 


627 6 



"From Ref. 22. 
''From Ref. 23. 



poor description of the ground state within the EXX ap- 
proximation. Using instead, our new potentials to recal- 
culate the mentioned properties results are in excellent 
agreement with experiment for all atoms studied. 

We can thus state with confidence that we now have 
at our disposal an affordable and not too complicated 
way of obtaining accurate results for a large number of 
static and low-frequency dynamic properties of physical 
systems. And all this within a relative minor modifica- 
tion of the EXX approximation within TDDFT. 

Of course, we have, at this stage, no idea about the pos- 
sible conserving properties of the proposed scheme. This 
is an interesting project for future research. But, on the 
other hand, conservation laws and sum rules might not 
be of such great importance to the properties discussed 
here. 

We should also mention that the results obtained here 
rely on the assumption that one can neglect the density 
and the frequency dependence of the XC kernel of the 
EXX. Based on the fact that this is no assumption in the 
case of a two-electron system we believe in the validity 
of this assumption but this must, of course, be further 
investigated. In the meanwhile wc propose to use the 
new method as an effective tool for calculating physical 
properties within the realm of TDDFT. 
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V. CONCLUSIONS AND OUTLOOK 



In the present work we have decided to temporarily 
part from our familiar, systematic and conserving way 
of constructing improved approximations to XC poten- 
tials and kernels within TDDFT in the linear regime. In- 
stead, we have tested the first step in an iterative scheme 
originally proposed by Peuckert. The starting point has 
been the previously studied EXX approximation for the 
XC potential and corresponding kernel. The method is 
described in detail in the sections above. The rational 
for this deviation has been the extraordinary accurate 
potentials obtained from this approach. The potentials 
are indeed very close to the exact XC potentials of DFT 
where these are known. Calculated ionization potentials 
for all atoms are close to experimental results in accor- 
dance with the well known fact that the highest occupied 
DF eigenvalue should equal minus the ionization poten- 
tial. 

We have in previous work calculated static polariz- 
abilities, low lying particle-hole excitation energies, and 
van der Waals coefficients from the EXX approximation. 
These turned out to be very similar to those of TDHF 
theory and not very accurate. We argued then that these 
less satisfactory results were a consequence of a rather 
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